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We show that with an appropriate surface modification, a slab of photonic crystal can be made 
to allow wave transmission within the photonic band gap. Furthermore, negative refraction and all- 
angle-negative-refraction (AANR) can be achieved by this surface modification in frequency windows 
that were not realized before in two-dimensional photonic crystals [C. Luo et al, Phys. Rev. B 65, 
201104 (2002)]. This approach to AANR leads to new applications in flat lens imaging. Previous 
flat lens using photonic crystals requires object-image distance u + v less than or equal to the lens 
thickness d, u+v ~ d. Our approach can be used to design flat lens with u+v — ad with a 2> 1, thus 
being able to image large and/or far away objects. Our results are confirmed by FDTD simulations. 

PACS numbers: 



I. INTRODUCTION 

Negative refraction (NR) was proposed theoretically a 
long time ago It was realized only recently in two 
classes of materials. One type is the so-called meta- 
materials @, [H, 0| , consisting of wires and split-ring res- 
onators. The second type is photonic crystals @, which 
have periodic permittivity and/or permeability. 

In these artificial materials, there are up to now two 
mechanisms to achieve negative refraction (NR). One is 
to use the anti-parallelism between the wave vector and 
the group velocity ((J. This can be realized in isotropic 
metamaterials |2| and in the second or even higher band 
around the T point in photonic crystals (PCs) [y]. In this 
case, the phase refractive index n p is negative. 

The second mechanism is to use the anisotropy and 
the concavity of the equi-frequency surface (EFS) such as 
the EFS around the M point in the first band of a square 
lattice PC. In this case though the phase refractive index 
rip is positive, the group refractive index n g is negative 
7]. In order to show negative lateral shift or all-angle- 
negative-refraction (AANR) , one has to orient the lattice 
such that the TM direction is along the surface normal. 
Negative refraction and flat lens imaging [8[ using both 
mechanisms have been observed in microwave ftl lloL [ill . 
[HI . ITiil , 0, EH and near infrared experiments [l6[ . 

Recently we proposed a new mechanism for negative 
refraction using surface grating [l7T ], This mechanism 
combines photonic band gap with surface grating to 
achieve NR and AANR. Negative lateral shift and AANR 
have been demonstrated in a multilayered structure with 
surface grating. 

In this paper, we obtain new windows of AANR in 
two-dimensional (2D) PCs using this new mechanism of 
surface modification. We also show that flat lens made of 
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photonic crystal with surface grating can have u + v = ad 

fwith it> 1 while for the Veselago-Pendry flat lens 
a = 1. Here d is the thickness of the flat lens, u and 
v are the distance from the lens to the object and the 
image, respectively. Thus a flat lens can focus large and 
far away objects. 



II. NEW APPROACH TO AANR 

In a pioneering paper [7| , Luo et al showed that within 
certain frequency windows in the first band of a PC, 
AANR can be achieved. Specifically, within the first 
band, AANR is possible only along the TM direction for 
a square lattice PC. We will show that with appropriate 
surface grating, NR and AANR is possible along the TX 
direction in the first band of a square lattice PC. 

In the main text of this paper, we only consider the 
transverse magnetic (TM) modes of a square lattice PC 
of square rods. Square lattice of circular rods or even 
rhombus rods can be treated similarly. The generaliza- 
tion to transverse electric (TE) modes and lattice struc- 
tures other than square lattice is straightforward. As a 
specific example, we consider a square lattice of square 
rods with size b/a = 0.7, thus a filling ratio of 0.49. The 
EFS is of this PC is calculated using the plane-wave ex- 
pansion [5| with 5041 plane waves. The EFS of the first 
band is shown in Fig. Q] The frequencies at the X 
point and the M point are ojx — 0.1943 x 2nc/a and 
u>M = 0.2446 x 2-Kc/a, respectively. 

Consider a slab of this PC with surface normal along 
the TX direction. If one increases the frequency, uj > u>x , 
there will be a partial band gap for waves incident to the 
air-PC interface since the Bloch waves have k y > k a and 
the incident plane wave with k y < k a will be completely 
reflected. Here k a is the k y value of the crossing point 
of the EFS with the XM boundary of the first Brillouin 
zone, as can be seen in Fig. Q] and Fig. [2j For certain 
frequencies u>i < u> < ujm with uji = 0.2089 x 2irc/a when 
k a = 0/c, a flat slab of such PC is an omnidirectional 
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FIG. 1: (Color online) EFS of the TM modes of a square 
lattice PC. The lattice is made of square rods of alumina 
(e = 9) with filling ratio 0.49 (see insert). The shading area is 
the window of AANR for the PC with surface grating a s = 2a 
(see Fig. [3j . A slab of this PC is oriented such that the surface 
normal is along the FX direction. 



reflector [lj|. For example for uj = 0.219 x 2nc/a as 
shown in Fig. [2l there will be total external reflection for 
any incident plane wave. However for these frequencies, 
a surface grating with period 
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(1) 



will give a momentum boost along the surface to the in- 
cident plane wave with incident angle 9 such that it will 
be coupled to the Bloch waves inside the PC with trans- 
verse momentum k y — ir/a + (0/c) sm9 if 9 is negative 
and k y — — ir/a+ (0/c) s'm9 if 9 is positive. The refracted 
wave will propagate on the opposite side of the surface 
normal with respect to the incident beam. Thus NR is 
achieved. This is illustrated in Fig. [2j The effect of this 
surface grating is equivalent to bringing down the EFS 
around the M point to the X point for u>i < u> < ojm- 
As we pointed out in Ref. [ItJ , it is the surface periodic- 
ity which determines the size of the EFS and the folding 
of the band structure. Furthermore, if ir/a — k a > 0/c, 
AANR can be achieved. The upper limit for AANR is 
lu u = 0.2192 x 27rc/a. Thus we obtained a 4.7% AANR 
around lu u . 

The above approach to NR and AANR is confirmed 
in our numeric simulations using FDTD [20]. Here we 
consider the lateral shift of an incident beam by a slab 
made of a square lattice PC. The detail of the slab is 
shown in Fig. [3] Negative lateral shifts are observed for 
different incident angles as shown in Fig. [4] for beams at 
oj = 0.219 x 2%c/a. It can be verified that for this slab, 
AANR can be achieved for 0.2089 < 0a/27rc < 0.2192. 
Note that the details of the surface grating is not essen- 
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FIG. 2: (Color online) Mechanism for NR and AANR for a 
square lattice PC with surface grating a s — 2a. The red curve 
is the EFS in the PC and the green one is that in the air for 
= 0.219 x 2-Kc/a. With the momentum kick Afc = n/a, an 
incident plane wave with fcj will be refracted negatively into 
a plane wave with k r . 



tial except its period a s — 2a. The grating can be holo- 
graphic or an array of circular rods as long as it is not 
too thick. For the specific surface grating shown in Fig. 
121 the energy transmissions are 99.2% and 4.7% for plane 
waves with incident angles 15° and 30°, respectively. 
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FIG. 3: (Color online) Details of slab made of a square lattice 
PC with surface grating a a — 2a. For simplicity, the thickness 
d of the slab is defined as the distance from the first surface 
to the last surface of the structure. 



III. FLAT LENS WITH LARGE a 

One prominent application of negative refraction is the 
Veselago-Pendry perfect lens [8j . A flat slab of thickness 
d can focus an object with distance u on one side to a 
distance v on the other side with u+v = d if the refractive 
index n = —1. For a generalized flat lens without optical 
axis [18], the lens equation takes the form 
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FIG. 4: (Color online) Negative lateral shift by a slab of PC 
given in Fig. [3] for incident Gaussian beam with incident 
angles 15° (top) and 30° (bottom) at = 0.219 x 2ixc/a. The 
distance is measured in the unit of lattice spacing a. 



with a a material property, depending on the dispersion 
characteristics of the flat lens. This is illustrated in Fig. 
[5j This lens equation requires the following form of the 
EFS at the operating frequency 



K — T^J UJ 2 /c 2 



V2 
K y . 



(3) 



The lens surface is in the y-direction. The surface normal 
is along the x-axis. Here k rx is the longitudinal compo- 
nent of the wave vector in the lens medium and k is the 
center of the EFS ellipse [ll]. 

Even though AANR can be realized in the first band 
along the TM direction for a square lattice PC [7|, the 
EFS is very flat around the lens normal. Thus one has 
er <C 1 Pi|. Though a ~ 1 has been reported in PCs 
with other structures [U, [22j , the focusing is still limited 
to the vicinity of the lens surface [23| . For practical ap- 
plications, we need large a so that the object and image 
can be far away from the lens. We will show that with 
our new mechanism for AANR, large a can be achieved. 




FIG. 5: A flat lens with lens equation u + v — ad. 



As we have stated in the previous section, the first 
Brillouin zone of a square lattice PC with surface grating 
a s = 2a takes the shape of a rectangle instead of a square 
and its vertical size is reduced to —n/a < k y < n/a. The 
center of the EFS for uii < uj < u>m of the original PC is 
moved from the M point to the X point. The fitting of 
the modified EFS for these frequency by Eq. ((3]) will give 
the lens property a. An inspection of the band structure 
reveals that the EFS is not elliptical. This results in 
incident angle dependent a [18|. Nevertheless, the EFS 
can be can be fitted well with a constant ao for small k y 
as shown in Fig. [51 For the square lattice PC we have 
designed (Fig. [3]), one has ao ~ 4 for u>i < uj < u u . 
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FIG. 6: (Color online) Fitting of EFS of the TM modes at 
0a/2ixc = 0.219 by Eq. ©. Here a = 4 and na/2nc = 1.16. 
Note that the center of EFS is shifted from the M point to 
the X point due to the surface grating a s = 2a. 

Focusing by such a flat lens is shown in Fig. [7) For 
a point source with u — 13.6a, a clear focused image is 
obtained at v = 12.4a for the operating frequency = 
0.219 x 2irc/a, which is consistent with the lens equation 
u + v = CTeffd with CT c ff = 3.5 and d = 7.4a. There are two 
reasons for cr ff < cr . First that the EFS is elliptical only 
for small k y = (0/c)sin# and a = ~dk rx /dk x decreases 
with increasing incident angle 9. The effective a e s is an 
average and thus smaller than ctq. Second, the thickness 
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of a PC slab is not a well-defined quantity. Here we 
simply define the lens thickness as the distance from the 
first surface to the last surface as shown in Fig. [3] This 
may overestimate the effective thickness of the lens. 
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FIG. 7: (Color online) FDTD simulation of flat lens focusing 
of a point source. For better contrast effect, the field intensity 
at the point source is suppressed. The detail of the lens is 
given in Fig. [3] 

To further check the performance of this flat lens, we 
vary the object distance. In Fig. [5] we show that the 
ratio (u + v)/d is almost constant and very close to do 
for different object distance u. 
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FIG. 8: The ratio of the object-image distance to the slab 
thickness (u + v)/d vs. the object distance u for the flat lens 
shown in Fig. |3]at the operating frequency = 0.219 x 2nc/a. 

The primary concern of practical application is the 
power transmission through the lens. As expected, the 
transmission through the flat lens is low due to the 
impedance mismatch. However since the details of the 
grating on the PC will not alter the scenario for NR, the 
power transmission can be enhanced through careful en- 
gineering of the grating. In our simulation, we find that 
the grating on the 2D PC with large dielectric constant 



has strong power transmission. The parameters of the 
surface grating shown in Fig. [3j are not optimized. Fur- 
ther improvement of transmission may be possible. For 
the grating parameters given in Fig. [31 the transmission 
coefficient is calculated and plotted in Fig. [9] 
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FIG. 9: (Color online) The transmission coefficient for plane 
wave incident on the flat lens shown in Fig. [3] 



IV. DISCUSSION AND CONCLUSION 

In this paper, we have achieved NR using a new ap- 
proach: photonic band gap with surface grating. This 
approach to NR gives new window for AANR in 2D PCs. 
This approach also gives flat lens made of these PCs with 
a large object-image distance. Thus these flat lenses are 
able to image large and far away objects. 

Our new approach gives much larger window of AANR 
than previous realized. For example for a square lattice 
air holes in e = 12 with r/a = 0.35 studied in Ref. [3], 
our approach gives a lower limit uji = 0.183 x 2nc/a and 
an upper limit lu u — 0.206 x 27rc/a, hence a fraction of 
AANR frequency range of 11% around 0.206 x 2wc/a. 
This range is much larger than the AANR range along the 
TM direction for the TE modes. This window of AANR 
is easier to locate. The two limits are obtained from the 
crossing of the band with the light lines around X point 
and M point, respectively. For the determination of the 
lower limit 0; , there is no need to compute the frequency 
at which the radius of curvature of the contours along 
TM diverges Q. 

Our new approach to AANR can also be extended to 
three-dimensional photonic crystals. 
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